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Abstract. We describe a simple method to study the dependence of the solar properties on a 
generic (small) modification the physical inputs adopted in standard solar models calculations. 



1. Introduction 

A new solar problem has emerged recently. The determinations of the photospheric heavy 
element abundances \1\ indicate that the sun metallicity is lower than previously assumed 
[2]. Solar models that incorporate these lower abundances are no more able to reproduce the 
helioseismic results. Detailed studies have been done to resolve this controversy (see [3] for a 
review), but the origin of the "solar composition problem" has still to be identified. 

In this situation, it is useful to perform a detailed analysis of the role of physical inputs and of 
the standard assumptions for standard solar model (SSM) calculations. This task is not always 
possible in simple terms, since the construction of SSM relys on (time-consuming) numerical 
integration of a non linear system of partial differential equations. In order to overcome this 
problem, we provide a tool which is at same time simple and accurate enough to describe the 
effects of a generic (small) modification of the physical inputs on the solar properties . The 
starting point is the idea that, despite the present disagreement with helioseismic data, the SSM 
is a rather good approximation of the real sun. We can thus assume small variations of the 
physical and chemical properties of the sun with respect to the SSM predictions and use a linear 
theory to relate them to the properties of the solar plasma. 

2. Linear Solar Models 

We consider the effects of a generic (small) modification of the input parameters of the sun. 
Namely, we indicate with 5n{r) and 5e{r) the relative variations of radiative opacity n{r) and 
of the energy generation e(r) along the SSM profile, i.e.: 



where I = k, e and r indicates the distance from the center of the sun. The notation Q 
indicates, here and in the following, the SSM value for the generic quantity Q. As a result of 
this modification, we obtain a solar model which deviates from SSM. We indicate with 5h(r) 
the relative variation of the physical quantities h = I, m, T, P with respect to SSM predictions, 
i.e.: 



I(T(r),p(r),Y(r),Xi(r)) = I(T(r),p(r),Y(r),Xi(r)) [1 + 5I(r)] 
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Figure 1. The coefficients Jh(f) (left panel), a' h (r) (middle) and f3' h {r) (right). See [4] for details. 



In the radiative region of the sun (r < = O.73Oi?0), the functions 5h(r) are described with 
good accuracy by the solutions of the linear system of ordinary differential equations: 
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[jp 5P + jt5T - 5m + AYi ni + % 5e] 
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have been calculated 



The coefficients 7/1, and a' h and the scale heights lh 
in [3] and are shown in fig. 1 and fig. 2 (left panel). The parameters AYjnj and 5C represent 
the absolute variation of the initial helium abundance and the relative variation of pressure 
at the bottom of the convective enevelope and can be univocally determined by imposing the 
appropriate integration conditions. At the center of the sun (r = 0) we have: 



5m = 7p j0 5P + 7t,o ST + 7^0 AYj ni + % fi 5eo 
51 = PpnSPo + faoSTo + PynAY^ + p'cQSC + p'odeo 
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where the subscript "0" indicates that a given quantity is evaluated at r = 0. At the bottom of 
the convective envelope (r = i?b), we have instead: 
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where A' Y = 0.626 and A' c = 0.025 and m c = M conv / 'Mq = 0.0192 is the fraction of solar mass 
contained in the convective region. 

The above equations completely determine the physical and chemical properties of the 
"perturbed" solar models. They have been obtained in [3] by expanding to linear order the 
structure equations of the present sun close to the SSM solution. We, thus, refer to the models 
calculated by using this approach as Linear Solar Models (LSM). In the derivation of eqs.Q we 
also assumed that: 



(i) The relative variation of density 5p{r) can be described by: 

5p(r) = 5P(r) - 5T{r) - P Y AY(r) 



(4) 



where -Py(r) ~ — dln/i/dY ~ — 5/[8 — 5Y(r)] and \i represents the mean molecular weight. 
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Figure 2. The scale heights lh(r) (left panel) and the coefficients £ft(r) (right panel). See [4] for details. 



(ii) The helium produced by nuclear reactions, Y nuc (r), is roughly proportional to the energy 
generation coefficient in the present sun , i.e. Y nuc (r) oc e(r). In this hypothesis (and 
neglecting the variation of elemental diffusion efficiency in the radiative region), the absolute 
variation of helium AY(r) is given by: 

AY(r) = £y (r) AYni + £ p (r) 5P(r) + £ T (r) 5T(r) + £ e (r) 5e(r) (5) 

The coefficients £/i( r ) have been calculated in |1] and are shown in fig. 2 (right panel). 

(iii) The global effect of elemental diffusion (i.e. integrated over the sun history) on the solar 
surface abundances is roughly proportional to diffusion efficiency at the bottom of the 
convective envelope in the present sun. In this assumption, one has: 

AY b = 0.838 AY^ + 0.033 5C (6) 

where the subscript "b" indicates that a given quantity is evaluated at the bottom of the 
convective region. By taking into account that (Z^/X^) is observationally fixed, we can 
also calculate the relative variation of metals in the convective region: 

5Z h = -1.088 - 0.043 5C (7) 

and in the radiative region (where we neglect the variation of diffusion efficiency): 

SZ(r) ~ 5Z ini = -0.887 AYni - 0.164 5C (8) 

(iv) The relative variation 5R\, of the convective radius Ry, can be calculated "a-posteriori" . By 
applying the Swartzchild criterion to the solution of eqs.Q and by expanding to first order 
(see [1] for details), we obtain: 

5R h = 0.449 AYni -0.117 5C - 0.085 5n h (9) 

In ref. [1] the predictions of LSM for helioseismic observables and neutrinos fluxes have been 
compared with the results of full non-linear SSM calculations for several selected cases, showing 
that a very good agreement is achieved. 
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